o 



o 



Analytical Approximations of 

Whispering Gallery Modes in 

Anisotropic Ellipsoidal Resonators 

Marco Ornigotti 1 
Andrea Aiello 1 2 

l Max Planck Institute for the Science of Light, Giinther-Scharowsky-StraBe 1/Bau24, 91058 

Erlangen, Germany 



Institute for Optics, Information and Photonics, University of Erlangen-Niirnberg, 



Staudtstra&e 7/B2, 91058 Erlangen, Germany 
marco.ornigotti@mpl.mpg.de 



Abstract: Numerical evolutions of whispering gallery modes of both 

isotropic and anisotropic spheroidal resonators are presented and used to 
i— i build analytical approximations of these modes. Such approximations are 

carried out mainly to have the possibility to have a manageable analytic for- 

' j—j mulas for the eigenmodes and eigenfrequencies of anisotropic resonators. A 

Q_i qualitative analysis of ellipsoidal anisotropic modes in term of superposition 

of spherical modes is also presented. 
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1. I. Introduction 

Thanks to their extremely high quality-factor (up to 10 10 ), whispering gallery mode (WGM) 
resonators are nowadays very promising devices for applications in optoelectronics |1|, bio- 
sensing (21 0, nonlinear optics [4] and other fields of applied physics. The interested reader 
on possible applications of WGM resonators is addressed to Ref. [ 1 1 and references therein. 
Although the first analysis of such devices dates back to 1939 in a paper by Richtmyer f5j, 
this topic remained unconsidered for many years until recent times, when the need of compact 
optical devices with great performances, together with a technological boost that gave the pos- 
sibility to manufacture not only silica microspheres but WGM resonators of any material (such 
as crystalline resonators) and shape (spheroidal J6), toroidal Q and so on), had contributed to 
make WGMs to experience a new renaissance. In these years, the theory of WGMs in micro- 
spheres was well established both analytically and numerically and allowed precise calculation 
of eigenmodes, radiative losses, field distributions, et cetera [ 1 ]. Unfortunately, when approach- 
ing the study of anisotropic devices (more interesting from the point of view of the applications 
due to their very good performances), analytical solutions cannot be found for geometries dif- 
ferent from an ideal sphere [ 14] or cylinder: in this case, it is more convenient to solve the 
problem with the help of numerical codes. 

Sometimes, however, an analytical solution (even if not so accurate as the numerical solution) 
would be preferred either because it is easier to manage or because one has the possibility to 
directly catch the physics of the problem by directly looking at its form. But how to validate an 
analytical approximation of a solution when a real analytical solution does not exist? To answer 
to this question, we propose in this paper to use a numerical finite element solver (COMSOL 
Multiphysics [ 16 1) to firstly find an appropriate approximation of the problem in exam and then 
use that approximation to check the validity of an analytical model for the solution. 

Although this method is potentially applicable to a vast class of problems, in this paper 
we focus our attention to a specific problem: finding whispering gallery modes (WGMs) in 
isotropic and anisotropic dielectric ellipsoidal resonators. Even if there are some works about 
the approximations of the WGMs for isotropic non-spherical bodies that present analytical ap- 
proximations in terms of spherical-like functions fTTl[T3l or superposition of spheroidal modes 
ifTOl . to the best of our knowledge no papers have been written on analytical solution for the 
anisotropic case. We then intend to fill this gap by finding with COMSOL a suitable approxi- 
mation for WGMs in these resonators and then presenting a correspondent robust and accurate 
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Fig. 1. (color online) Sketch of the 3D geometry of the resonator (left) and its cross section 
in the xz plane (right). The numerical values of the semiaxes are a = 1.9 mm and b = 0.7 
mm.The cross section on the xy plane is a circle of radius a. 



analytical solution. The approximation we present is based on the observation that since WGMs 
are localized near the resonator surface, an effective approximation of an ellipsoidal resonator 
near its surface could be represented by a toroidal resonator of circular cross section, the radius 
of the torus being equal to the major axis of the resonator and the radius of the circular cross 
section being equal to the rim radius of the resonator at the considered surface. 

Within this framework we calculate both the field distributions and the eigenvalues of the 
WGMs of both ellipsoidal resonator and toroidal resonators, showing how the WGMs of the 
torus well approximate the WGMs of the ellipsoidal resonator. This result gives us the possi- 
bility to substitute, under certain conditions, the complete set of spheroidal wave functions ifTUl 
that characterize the ellipsoidal resonator with the complete set of spherical wave functions 
|[T4l[T5l . This permits us to manage a simple and analytical set of eigenmodes that can be easily 
used for theoretical predictions of nonlinear or quantum optical effects in these resonators. 

This paper is organized as follows: in section II the model used to simulate with COM- 
SOL such a resonator is presented, together with a brief recall of the weak-form expression for 
Maxwell's equations in a rotationally-symmetric resonator that is used as a model for the calcu- 
lations. Section III is then devoted to present the results of the COMSOL simulations we made 
both for the isotropic and anisotropic case, with a direct comparison between the field distribu- 
tions and the eigenvalues of both the ellipsoidal and toroidal resonators. Finally, in section IV 
conclusions are presented. 

2. II. The Model 

2.1. A. Geometry of the resonator 

Let us consider an ellipsoid of revolution filled by a dielectric medium, rotationally invariant 
around the z-axis and characterized by the in-plane major axis a, i.e. the radius of the circum- 
ference in the plane xy, and the out-of-plane minor axis b, i.e. the minor axis of the ellipse 
in the plane xz. Such an ellipsoid is depicted in FigJT] and its cartesian equation is the following: 
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where the major axis is a = 1.9 mm and the rim radius p = 0.25 mm; the value of the minor 
axis is then determined from these two numerical values with easy geometrical considerations 
and it turns out to be b = 0.7 mm[22 1. Following Ref.[4], we consider this resonator to be made 
of Lithium Niobate (LiNbOs), an anisotropic crystal characterized by the following dielectric 
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Fig. 2. (color online) Sketch of the 3D geometry of the toroidal (left) resonator that ap- 
proximates the ellipsoidal resonator and its cross section in the xz plane. The radius of the 
circular cross section is r = p = 0.25 mm, and the external radius of the torus is equal to 
a= 1.9 mm. 



tensor: 



£ = 




(2) 



where e± = 5.3 is the dielectric constant in the xy plane and £ii = 6.5 is the dielectric constant 
parallel to the z axis. However, in considering the resonator as isotropic, we assume that its 
dielectric tensor will be diagonal, with nonzero elements that are equal to e± . 

In order to make a suitable approximation, we compare the results obtained from this res- 
onator with the ones obtained from a toroidal resonator with circular section like the one de- 
picted in Fig. p] Here the external radius of the torus being a = 1 .9 mm equal to the major axis 
of the ellipsoidal resonator, and the radius of the circular cross section being r = p= 0.25 mm, 
equal to the rim radius of the ellipsoidal resonator. A sketch of this approximation is shown in 
Fig. [3] where both the lateral and top sections of the two resonators are shown. 

Even if an analytical solution for the open ellipsoid (with the fields continuous at the 
resonator surface) does not exist, it is possible to transform the open resonator into a closed one 
(with ideal metallic boundaries, i.e., the electric field is zero at the resonator surface) by simply 
considering that the field of a WGM outside the resonator is evanescent (i.e. exponentially 
decaying from the resonator surface) and replace the real spheroid with semiaxes a and b with 
effective semiaxes a = a + a and b = b + a, where a is the depth upon which the optical field 
penetrates in the surrounding medium under total internal reflection. For grazing angles, <7 
does not depend on the angle and has the value: 
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where ka is the vacuum wavenumber, % = 1 for quasi-TE modes and % = e for quasi-TM mode, 
being £ the dielectric constant of the isotropic resonator. A detailed analysis on the origin of 
this term can be found in Ref. [ 13 1. For the rest of the paper we will implicitly assume that the 
resonator is closed with semiaxes a and b, remembering that the solutions we will present can 
be easily adapted to the open resonator by simply substituting a and b with a and b. 
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Fig. 3. (Color Online) Lateral (a) and top (b) view of the ellipsoidal resonator (blue solid 
line) and the approximating torus of circular cross section (red shaded circle). The radius 
of the torus cross section r matches the rim radius p =b ja of the ellipsoid. The external 
radius of the torus matches the major axes of the ellipsoid (that appear as a circle of radius 
a in the top view), while the internal radius of the torus is such that its circular cross section 
has a radius of r. 



2.2. B. Weak formulation of the electromagnetic problem 

Before entering deep in the subject of this paper, we briefly intend to recall to the mind of the 
reader the weak formulation of Maxwell's equations with Galerkin's method of the weighted 
residuals, that we use to simulate with COMSOL whispering gallery modes in a general ax- 
isymmetric dielectric medium with permittivity tensor e. 

We choose to implement such a model rather than using a standard mode solver because the 
latter cannot be easily configured to fully exploit the axial symmetry of the problem and they 
experience some problems when dealing with WGMs. Then in order to obtain from them an 
accurate solution a fully 3D model must be implemented, making the calculations very time 
consuming and complicated. 

Galerkin's method, on the other hand, has the advantage that can easily take into account 
the symmetry of the problem, giving in this case the possibility to reduce the dimension of the 
problem from 3D to 2D, by solving the problem only in the xz plane thanks to the axial rotation 
symmetry along the z axis of the resonator. 

This method, compared with mode solvers, allow us to save a lot of computational time and 
memory, and gives us the possibility to calculate all the fields in the post processing phase of 
the simulation and only if we are interested in those quantities. 

This method can be applied either for solving directly the electric field components |20| or 
the magnetic field components ETI . depending on what is the best choice with respect to the 
problem to be solved. 

Our results are based on the discussion presented in Ref. [20|, where the calculation of the 
electromagnetic field distribution on an axisymmetric resonator is generalized to an arbitrary 
number of dielectric media embedded in an outer dielectric resonator, and they will be par- 
ticularized in this work to study the simpler case of a resonator made by only one dielectric 
medium. 

The electromagnetic field inside the resonator obeys Maxwell's equations in continuous 
macroscopic media [17). For this reason, in general, if one assumes that the resonator's con- 
stituent medium has constant magnetic susceptibility, then the magnetic field H will be continu- 
ous across the resonator volume, and the problem can be easily solved in terms of the magnetic 
field rather than the electric field E. 



Vx (r l VxH)-«v(v.H)+ 1 ^H = 0, (4) 



In the general case, however, the medium is not isotropic, and the dielectric constant be- 
came a tensor. For the sake of clarity, we will approach the problem from a general point of 
view, leaving the dielectric constant as a tensor in order to obtain a general expression for the 
magnetic field inside such a resonator. 

After obtaining the general expression, we will specialize it to either the case of isotropic 
(Section III. A) or anisotropic (Section III.B) resonator. 

By eliminating the electric field in favor of the magnetic field in the Maxwell's equations, 
after some simple algebra the problem reduces to solve the following equation: 

]_d 2 
c 2 

where a is the so-called penalty constant that controls the presence of spurious solutions (i.e. 
solutions with nonzero divergence inside the resonator) that may arise during the computation. 
The presence of this term is quite common when a weak form differential problem has to be 
solved with the help of finite element codes (see for instance Ref. |18|). For all the simula- 
tions presented in this paper, the penalty constant is equal to one. £ _1 is the inverse relative 
permittivity tensor (sometimes named in literature as the impermeability tensor), assumed to 
be independent of field strength. 

In order to have the complete solution of this equation, suitable boundary conditions must be 
applied: in our case we will assume that the resonator is closed and so we will use the so-called 
electric wall boundary conditions, namely: 

H n = 0, (5a) 

E x n = 0. (5b) 

These two equations imply that the magnetic field must be tangent to the resonator surface, 
while the electric field is normal to the resonator surface. 

From Eq. Q, and with the aid of Galerkin's method of weighted residuals 1 19 1 it is possible 
to obtain the desired weak form of the electromagnetic problem, that will be implemented in 
COMSOL to be solved. 

In order to reduce Eq. Q in its weak form counterpart, we multiply both sides of Eq. Q 
by the complex conjugate of a test magnetic field H* and then integrate over the resonator's 
dielectric volume. 

Then we integrate by parts over the spatial coordinates and after the disposal of the surface 
integrals with the help of boundary conditions, we finally arrive to the desired weak form of 
Eq.Q, that reads: 



/[(vxfl*).(IvxH)-a(v.H*)(v.H^H 
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dt 2 



dV = 0, (6) 



where the integral is taken over the whole resonator volume, and the first term under the integral 
stands for £? =1 [V x H*],e ( ^ [V x H] 7 . The three terms appearing in the integrand correspond 
exactly to the weak-form terms required to define an appropriate model within a partial differ- 
ential equation solver [ 19 j. 

The axial symmetry of the problem suggests to describe the resonator in cylindrical coor- 
dinates {r,z,(p}- With this choice of coordinate system the above equation turns into a 2D 
equation, because the ^-dependent part of the magnetic field H can be factored out in the 
form of a complex exponential exp(iM(p), where M is the azimuthal quantum number, i.e.. 
the number of nodes of the WGM in the xy plane. The total magnetic field can be then writ- 



ten as H(r) = exp(/M<p) H r (r,z)f + H z (r,z)z + iH<p(r,z)(f> 



where the imaginary unit in the 
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Fig. 4. (color online) Contour plot of the fundamental isotropic WGM (m = 1 000, / — m = 0, 
q = 1) of (a) the ellipsoidal resonator and (b) its toroidal approximation. The dashed black 
lines represent the resonators boundary. 



^-component is added to allow writing all the three components of the magnetic field in terms 
of an amplitude multiplied by a common phase factor. As it can be seen, the initially full 3D 
problem has been reduced to a 2D problem, thus resulting in a simplification of the problem 
and a speeding up of the computational time needed to execute the calculations. 

The next step is to transform Eq. |6]) and its solution to a form easy to implement in COM- 
SOL; this procedure is clearly presented in detail in Ref.|20|, so we address the interested 
reader to this paper. We just point out that for the isotropic case one has to put e± = £11 in Eq. 
(8) and following in Ref. [20] , being e± the dielectric constant in the xy plane of the resonator 
and £ii the dielectric constant parallel to the z-axis. In the case of uniaxial resonator, conversely, 
these two quantities have to be different. 



3. III. Results 

3.1. A. Isotropic resonator 

In Figs. |4]and|5]two examples of the field distribution for WGMs in the cross sectional xz plane 
for both spheroidal (left) and toroidal (right) resonator are shown. In particular, Fig. [4] shows 
the fundamental WGM characterized by the three quantum numbers I = m = 1000 and q = 1, 
being /, m and q the orbital, azimuthal and radial quantum numbers respectively. FigBlinstead 
shows the sixth WGM characterized by having I — m = 4 and q = 2. 

As can be seen from these figures, and from their x- and z-sections shown in Figs. 6]|9 the 
toroidal approximates the isotropic spheroidal resonator very well, resulting in the possibility 
to describe the set of spheroidal modes with the equivalent set of toroidal modes. Since here 
we are analyzing a 2D problem, and all that matters are the shapes of the modes in the cross 
section plane xz, saying that the toroidal modes well approximate the spheroidal mode is equal 
to say that these modes are very well approximated by the eigensolutions of the torus. 

In order to make the approximation complete, we calculated the eigenvalues corresponding 



to the resonator's eigenmodes. Fig. 10 shows a comparison between the eigenvalues of the first 
few WGMs of both spheroidal (red) and toroidal (blue) resonators. As can be seen, also in this 
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Fig. 5. (color online) Contour plot of the sixth isotropic WGM {m = 1000, l — m = 4,q = 2) 
of (a) the ellipsoidal resonator and (b) its toroidal approximation. The dashed black lines 
represent the resonators boundary. 




Fig. 6. (color online) Comparison of the x-section (radial part of the WGM) of both el- 
lipsoidal (red curve) and toroidal (blue curve) fundamental WGM. The black dashed line 
represents the resonator boundary. 



case the approximation holds very well. For the analysis and comparison of the eigenvalues of 
these two cavities, we adopted as a merit parameter for this approximation to hold the mean 
error on the estimation of the spheroidal eigenvalues with the toroidal ones (i.e. the numerical 
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Fig. 7. (color online) Same as Fig.[6]but for the sixth WGM. 
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average of the error on every single eigenvalue). This mean error takes a value roughly around 
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Fig. 9. (color online) Same as Fig.IHlbut for the sixth isotropic WGM. 
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0.3% for the isotropic resonators, i.e. on the order of 0(l~ l ) for the WGMs presented here. 
Normally one can find in literature [23 1 that the most common level of accuracy in determining 
the eigenvalues of a resonator is of 0(l~ l ) for WGMs. In our case this means that the precision 
needed for the eigenvalues must be of the order of 0.1%, that is precisely the accuracy lever of 
the toroidal approximation. 

The approximation is very good for the first few WGMs, where the mode is highly confined 
near the resonator surface both in x and z direction. This results in a complete equivalence 
between the ellipsoidal and toroidal boundary condition. As can be seen from Fig. 10 the 



eigenfrequencies are very close to each other, resembling the fact that the modes have a very 
good overlap. This can be justified by also saying that while the modes of both resonators are 
highly confined near the resonator surface, they experience the same boundary conditions, i.e. 
they are almost identical. 

When higher order modes are considered (like for example the one depicted in Fig.|9]l, some 
discrepancies start to appear, since the mode starts to penetrate deep in the resonator. Formally, 
this approximation breaks when the considered mode is too extended in the resonator (namely 
along the ^-direction) in such a way that its confining region is greater with respect to the region 
in which there is a good overlap between the curvature of the ellipse and the curvature of the 
circular cross section of the torus. In this case, the modes of the two cavities feel a different 
boundary (one feels ellipsoidal boundaries while the other one feels circular boundaries), and 
the discrepancy between these two modes starts to be significative. 

3.2. B. Anisotropic resonator 

For the case of anisotropic resonators, results are reported in FigfTTIand 12 for the bidimen- 
sional xz cross sections of both spheroidal and toroidal modes and in Figs 13fT6 for their sec- 
tions along the x- and z-axes respectively. 
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Fig. 10. (color online) Comparison between isotropic ellipsoidal (red dots) and toroidal 
(blue dots) eigenvalues corresponding to the first seven WGM. The inset shows the expres- 
sion of the relative error used in this work to estimate the validity of our approximation. 
Xt is the eigenfrequency of the torus, while Xe is the eigenfrequency of the ellipsoidal 
resonator. 



The geometry of the anisotropic resonators is the same as the one described before for the 
isotropic case. The value of the dielectric function in the xy plane and the one parallel to the z 
direction are assumed to be e± = 5.3 and £11 = 6.5 respectively, as stated in section II.A. 



Figure 17 shows instead the eigenvalue comparison for the anisotropic case for the first seven 
WGMs. As can be seen, even in the anisotropic case the toroidal approximation very well 
reproduces the eigenvalue pattern of the spheroidal resonator, and in this case the mean error 
was estimated to be roughly around 0.4%. 

3.3. C. Analytic formulas 

The approximations shown in the previous sections are of great interest, since they give us the 
possibility to approximate with a very good level of accuracy the WGM solutions of Helmholtz 
equations for an ellipsoid, by means of the WGM solutions of the Helmholtz equation in a 
toroidal cavity with circular cross section. This is very important from the theoretical point of 
view, because, as the WGM eigenmodes of the torus approximate well the WGM eigenmodes 
of the ellipsoid, then it is possible to use simple analytic expressions (as, for example, the 
solutions of the Helmholtz equation for the circle) to describe both isotropic and anisotropic 
ellipsoidal resonators. 

The results presented here show that if we take a circumference whose radius is equal to 
the rim radius of the actual spheroid near its boundary, the solutions for this geometry very 
well approximate the xz cross sections of the spheroidal solutions. In order to describe also the 
rotational symmetry of the spheroid, this circumference must be rotated around the z-axis, gen- 
erating a toroidal resonator like the one depicted in Fig. [2] For this resonator, the ^-dependent 



term in the eigenmodes simply consists (due to the axial symmetry) in the complex exponential 
of the form exp (iM(p), where M is the number of nodes of the eigenmode in the xy-plane, i.e. 
around the z-axis. This ansatz highly simplifies the analytical managing of such modes in these 
resonators, giving the possibility to work with the complete set of circular solutions (Bessel 
functions along the x-direction and trigonometric functions along the z-direction) rather than 
the more complicated spheroidal wavefunctions, concurring in a huge simplification of analyt- 
ical calculations in such a system. Such calculations usually represent very hard tasks, as they 
require to manage infinite series of spherical solutions |[T0l . It has to be pointed out, however, 
that this kind of approximation is not really needed when dealing with isotropic resonators, 
since a variety of analytical formulas is already available ifTTl [T21 [131 . The present approxima- 
tion, on the other hand, is of great importance for anisotropic resonators, for which an analytical 
expression for the eigenmodes does not exist. 

Having this goal in mind, let us first consider the radial wave equation for the circle, that has 
the following form: 

d 2 \\l 1 d\\l f,l 2 \ 

dx 2 x dx V x 2 J 
where x = kr. The solutions to this equation are the Bessel functions of the first kind 7;(x), 
as they are the only ones that are finite at the origin. Nevertheless, as usual, it is possible to 
approximate the Bessel functions of large argument with the correspondent Airy functions 
Ai(x) fl5|. By exploiting these large-limit solutions, it is possible to build an analytical 
approximation to the ellipsoidal eigenmodes in the form of suitable combinations of Airy 
function as follows: 



,v 



¥(x) = £ CfcAi Y[x-P+x k 

lc=0 



(8) 



where Xk is the k-th zero of the Airy function and the quantities 7 and j3 are two scaling and 
fitting parameters, respectively, that are used to take into account both the effects of anisotropy 
(7 parameter) and other displacement effects (/3 parameter). The weighting coefficients c\ has 
to be chosen in such a way that the sum can be faithfully truncated to the order N. With this 
expansion, it is then possible to address the reals WGMs of the ellipsoidal resonator with a 
finite number of solutions of the circle with large argument. 

With these expressions for the WGMs it is then possible to fit the ellipsoidal modes. For 
example, the fundamental WGM described in Fig.[l3]can be very well reproduced by using only 
one term in the expansion ([8]) with 7 = 25 and j3 = 4.62. In this case the term x^ — —2.338 1 will 
be the first zero of the Airy function of the first kind Ai(jt). For the more complicated pattern of 
Fig.[l4j instead, the ellipsoidal WGM can be very well approximated by using three differently 
weighted Airy functions, so that N = 3, with 7 = 23 and j3 = 4.46. The weighting coefficient 
in this case are c\ = 0.15, c^ = —0.62 and C3 = 1.8. 

Proceeding along this way, we have then the possibility to approximate, with a large degree 
of precision, all the WGMs of the ellipsoidal resonator (either isotropic or anisotropic) with a 
finite number of eigensolutions of the Helmholtz equation in the circle, rather than using infinite 
series. 

This is of great advantage because it guarantees more manageable expressions, and allows 
one to easily understand the physics that lives behind these modes, as one can interpret them as 
just a linear combination of few modes from a circle. 

In the same manner it is possible to analytically approximate the angular distribution of such 
modes, i.e., the WGM structure along the z direction. For this case we have the possibility to 
employ the solution to the angular equation for the circle 



d 2 © 



P- 



r® = o, (9) 



dz 2 
by writing again the angular part of the mode of the ellipsoid as a finite sum of trigonometric 



terms: 



N r 



®(z) = E b k cos(pkz)+d k sm(pkz) 

k=\ 



(10) 



where /3 = la is the fitting parameter (I being the orbital quantum number that appearing in the 
angular equation) that, again, takes into account a possible anisotropy even along the angular 
distribution, and b k and d\ are weighting coefficients. 

It has to be pointed out, however, that while for an uniaxial anisotropy in a sphere the 
anisotropy will affect (at the first order) just the radial part of the mode [14|, in the case of 
an ellipsoid we must account for an anisotropy fitting parameter even in the angular part of 
the mode since, as the mode index grows, the mode starts to extend itself into the resonator, 
breaking the toroidal approximation, then starting to feel ellipsoidal boundary conditions rather 
than the circular ones. This can be considered and modeled as an effective anisotropy that acts 
prevalently along the angular direction. 

For example, for the sixth mode of the anisotropic ellipsoidal resonator as the one shown 



in Fig. 16 for which / = 1000 and I — m = 4, a very good analytical approximation can be 



built with the help of Eq. ( fTO) by considering only the even part of the expansion (i.e., only 



cosine terms), because the mode shown in Fig. 16 possesses even symmetry with respect to 
the vertical axis. We then have N = 4 and the expansion coefficients are b\ = 2.6, hi = —1.5, 
£>3 = 0.1, bi, = —0.2; the effective anisotropy parameter takes the value a = 0.022 and d^ = 
since we consider only angular function with even parity. 

As can be seen, the value of the effective anisotropy parameter in the argument of the an- 
gular functions is very low but still different from zero. This means that while for a spherical 
anisotropic resonator (if the anisotropy is small) the angular function is not affected by the 
anisotropy, in the case of the anisotropic ellipsoidal resonator, this is also true for the low order 
WGMs. As the order of the WGM grows, than an effective anisotropy (arising from the dif- 
ferent shapes of the spherical and ellipsoidal boundaries that the mode sees), starts to play an 
important role even in the angular distribution of the mode. 

4. IV. Conclusions 

In this work we have used a numerical finite element solver like COMSOL to find the field dis- 
tribution of WGMs in a spheroidal resonator. We have proposed to approximate the spheroidal 
resonator near its boundaries with a toroidal resonator with circular cross section, having the 
radius equal to the major axis of the spheroid and the radius of the circular cross section equal 
to the rim radius of the spheroid neat its surface. We presented results for both isotropic and 
anisotropic resonators and we pointed out that, especially for anisotropic resonators, thanks 
to the numerical verification of our approximation, it is possible to approximate WGMs of a 
spheroid with a suitable superposition of modes of the circle. The advantage of our proposed 
analytical approximation is two-fold: from one side it gives the possibility to easily represent 
the ellipsoidal WGMs as superposition of few circular eigenmodes rather than infinite series 
of spherical harmonics. On the other side, for the case of anisotropic resonator, it gives easy 
manageable formulas to describe modes that otherwise will not have an analytical representa- 
tion. This could be of great help when challenging complicated calculations on such resonators 
because at first it allows us to use analytical (i.e., simple and practical) methods to have a first 
educated guess about how things should go inside the considered system. Second, a posteriori, 
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Fig. 11. (color online) Same as Figplbut for the anisotropic case. 
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Fig. 12. (color online) Same as Fig.Blbut for the anisotropic case. 




availability of simple approximate formulas, permits very powerful qualitative analysis of the 
quantitative numerical results. 
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Fig. 13. (color online) Same as Fig.[6]but for the anisotropic case. 
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Fig. 14. (color online)Same as Fig.[7]but for the anisotropic case. 
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Fig. 15. (color online) Same as Fig.lSlbut for the anisotropic case. 
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Fig. 16. (color online) Same as Fig.RJlbut for the anisotropic case. 
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Fig. 17. (color online) Comparison between anisotropic ellipsoidal (red dots) and toroidal 
(blue dots) eigenvalues corresponding to the first seven WGMs. The inset shows the expres- 
sion of the relative error used in this work to estimate the validity of our approximation. 
Xj is the eigenfrequency of the torus, while Xe is the eigenfrequency of the ellipsoidal 
resonator. 



